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INTRODUCTION 


Part V of this Final Report on "Improved Fluid Dymanics 
Similarity, Analysis and Verification" is entitled ANALYTICAL 
AND EXPERIMENTAL STUDIES OF THERMAL STRATIFICATION PHENOMENA. 

It consists of three sections, based on a Master's thesis 
and two publications respectively . Because of this, each of 
the sections is self contained and has its own page numbering 
system. The titles of these sections are: 

1. "Transient Stratification of Fluid in a Container 
Due to a Sudden Increase in Wall Temperature" 

2. "An Interferometric Study of Thermal Stratification 
by Natural Convection in a Contained Liquid" 

3. "Use of Multiple Outlets for Reducing the Stratification- 
Induced Propellant Supply Temperature Rise" 

In the first section an analytical study of transient 
stratification in a rectangular container is performed using an 
integral method for the case in which two of the vertical walls 
are heated and are maintained at a constant temperature, with 
the remaining container walls being insulated. Experimental 
measurements obtained with water as the test fluid are presented. 
These data compare favorably with the analytical predictions for 
much of the transient response period. Sizeable deviations 
occur during the latter portion of the transient period, however, 
due to limitations of the analytical model in accurately describing 
the actual physical phenomena in this regime. 



In the second section an experimental interferometric study 
of thermal stratification is described. In this investigation 
interferometer photographs were obtained. From these photographs 
the entire instantaneous temperature distribution can be 
visualized for any time during the transient stratification period. 
Temperature profiles obtained in this manner are compared with 
analytical predictions in this section/ and good agreement between 
the two is illustrated. 

In the third section a new technique is proposed for reducing 
the increase in the propellant supply temperature risO due 
to thermal stratification in the supply tank. The technique 
is based on proper control of the flow distribution through 
multiple outlets. By reducing the supply temperature rise which 
typically occurs during the latter portion of the engine burning 
period, the potential problem of cavitation in the supply pumps 
during this period is alleviated. Experimental measurements are 
given to illustrate the value of the proposed method. 



SECTION 1 


TRANSIENT STRATIFICATION OF FLUID IN A CONTAINER 


DUE TO A SUDDEN INCREASE IN WALL TEMPERATURE' 


^This section of Part V of the Final Report consists of a Master's 
thesis entitled "Transient Stratification of Fluid in a Container 
Due to a Sudden Increase in Wall Temperature",, by Terry Lynn 
Libby, School of Mechanical Engineering, Purdue University, 
Lafayette, Indiana, 
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ABSTRACT 


This investigation was undertaken to determine the stratification 
behavior of a contained fluid subjected to transient free convection 
heat transfer. The motivation for this study was provided by the fact 
that stratification processes are encountered in cryogenic fluids used 
in liquid propellant rockets. . For convenience, water was u$ed as the 
test liquid in conducting the experiments. A rectangular vessel was 
employed with heat transfer from two opposite walls of the vessel to the 
fluid. The wall temperature was increased suddenly to initiate the 
process and was then maintained constant throughout the transient 
stratification period. Thermocouples were positioned on a post at the 
center of the vessel. They were adjusted so that temperatures could be 
measured at the fluid surface and at specific depths beneath the 
surface . 

A major goal of this work was the attainment of analytical and 
experimental results for the surface temperature and for the stratified 
layer thickness. The predicted values of the surface temperature and 
the stratified layer thickness were found to agree reasonably well with 
the experimental measurements. The experiments also provided information 
on the transient centerline temperature distribution and the transient 


flow distribution. 
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INTRODUCTION 


When a fluid contained in a vessel is subjected to external heat- 
ing, a process known as stratification takes place. In this process 
the warmer fluid tends to rise to the top of the vessel due to its 
buoyancy, and this causes the fluid to be formed into horizontal layers 
having different temperatures. For this particular study the heat 
transfer occurred through two opposite walls of a rectangular vessel, 
as shown in Figure 1. 

The heat transfer through the container walls results in the 
creation of a free convection boundary layer along each wall. As 
the fluid in the boundary layer becomes heated it rises and is simulta- 
neously replaced by the colder bulk fluid near the bottom of the vessel. 

When the fluid rising along the wall reaches the surface, it leaves 
the boundary layer and flows across the surface toward the geometric 
center of the vessel. The flows from both boundary layers meet at 
the center of the symmetrical vessel. Here they form a plume which 
moves downward into the bulk liquid region. This downward motion of 
the plume is arrested rapidly and is reversed by buoyancy effects. 

The warmer fluid in the plume then moves upward and outward simulta- 
neously, forming a nearly uniform layer of fluid called the stratified 
layer. As the process continues the surface temperature increases 
with an ever decreasing rate and approaches a maximum value. Simultaneously, 
the stratified layer thickness approaches the total fluid depth with 
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a decreasing rate. A more complete explanation of the stratification 
process is given in the following section of this report. For purposes 
of the present discussion the brief description given above will suf- 
fice. 

The stratification of cryogenic fluids is very important in situa- 
tions where these fluids are used as rocket propellants. In these 
cases the bulk liquid temperature is very low, and large temperature 
differences typically exist between the bulk liquid and the environment. 
This causes large heat transfer rates to the fluid which can cause a 
sizeable fraction of the fluid to become heated and stratified. When 
this heated fluid is pumped rapidly to the engine during flight, its 
high temperature can cause cavitation to occur. A condition such as 
this in the pumps can choke off the propellant flow to the engines , 
and thus lead to premature engine burnout. 

The significance of stratification is not limited to the space 
industry. The stratification of a lake is beneficially put to use by 
knowledgeable fishermen. It is a verified fact that different fish 
desire different water temperatures. These temperatures are created 
by the stratification of the lake's water, and a knowledge of this 
phenomenon is a helpful aid to the fisherman. 

In the present investigation the fluid was contained in a rectan- 
gular vessel which permitted heat transfer through only two parallel 
walls (Figure 1) . The remaining walls were insulated. The stratifica- 
tion process was experimentally studied, and the measurements were 
compared with the corresponding predictions obtained by analysis. For 
convenience, water was used as the test liquid. A constant wall tem- 
perature was imposed at the heated surfaces. Experimentally, this was 
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obtained by suddenly flowing hot water at a preset temperature over the 
outside of the two walls of the vessel. 

An approximate integral method was used in the analysis to obtain 
equations which govern the stratification process. The use of this 
method implies that solutions of the resulting equations are not exact 
in their description of the stratification process. The analytical 
results are therefore not quantitatively precise, but they are qualita- 
tively correct. For this reason they serve very well the purpose of 
providing good approximate results of the process. 

In the following section the stratification process is described 
in further detail. Then a review is given of previous articles and 
reports describing investigations which have been conducted on the 
phenomena of stratification in order to provide sufficient background 
material. The specific problem investigated in this study is treated 
in the remainder of this presentation. 
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DETAILED DESCRIPTION OF THE STRATIFICATION PROCESS 


In this section a detailed physical description of the stratifica- 
tion process is given. The various features which are discussed here 
are well known from previous investigations . Some of these previous 
studies are reviewed in the following section of this report. 

As stated in the Introduction, the heat transfer through the two 
walls creates a free convection boundary layer along each wall adjacent 
to the body of stationary fluid. The fluid in these boundary layers 
becomes heated and consequently decreases in density. This decrease 
in density, together with the gravitational field, creates an upward 
buoyancy force and flow within the boundary layers. The fluid moving 
up and finally out of the boundary layers is supplemented by fluid 
moving into the boundary layers in the bottom portion of the vessel. 

This latter fluid which continuously replenishes the boundary layers 
comes from the bulk fluid. The heated fluid continues up the entire 
height of the boundary layer to the surface. This hot fluid then con- 
tinues across the surface in a direction normal to the heated walls 
and inward toward the center of the vessel. The flows from the two 
boundary layers meet just below the liquid surface at the center of 
the vessel. At this point they form a plume which moves downward toward 
the bulk fluid. The initial phase of the stratification process is 
depicted in Figure 2a. The process continues as the downward motion 
of the plume is arrested and is finally reversed due to buoyancy effects. 
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The warmer fluid within the plume then moves upward and outward simul- 
taneously, forming a layer of fluid of nearly uniform thickness called 
the stratified layer. At the same time, a reverse shear layer is formed 
as shown by the darkened area impinging on the boundary layer (Figure 2b) . 
The vorticity within this layer tends to induce downward motion in the 
rising fluid as well as motion across the fluid bulk under the strati- 
fied layer. As the process proceeds, the reverse shear layer becomes 
less pronounced due to its gradual suppression by the growing stratified 
layer. Likewise, the part of the boundary layer congruent with the 
stratified layer becomes very thin, as shown in Figure 2c. In this 
manner the entire fluid bulk eventually becomes stratified. As long 
as the boundary layers are maintained, the process continues until the 
entire fluid bulk reaches a uniform temperature equal to the temperature 
of the container walls. 

As a final note, the bulk fluid receives some heat by conduction 
from the stratified layer. An order of magnitude analysis shows that 
this effect is very small for time periods which are typical of the 
stratification process. 
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PREVIOUS INVESTIGATIONS OF STRATIFICATION PHENOMENA 

The stratification process has been studied by a number of inves- 
tigators. An experimental investigation of the stratification of a 
cryogenic fluid contained in a vessel was carried out by Scott et al. 
(Reference 1) . It was found that the presence of a good heat conducting 
rod placed vertically within the liquid provided a path for heat con- 
duction and eliminated the stratification tendencies of liquid helium 
contained in a nonventing Dewar vessel. 

Anderson et al. (Reference 2) have given an approximate mathemati- 
cal method for the analysis of the stratification of a contained fluid 
for the situation in which the vessel walls are heated by radiation. 
Barakat (Reference 3) , on the other hand, presents an exact mathematical 
method for the description of the flow patterns during the stratifica- 
tion process. 

A literature survey, carried out as a portion of the present study, 
revealed that most of the investigations conducted on stratification 
phenomena employed very specific conditions. Primarily, these condi- 
tions either involved the method by which the vessel walls were heated 
or were concerned with the geometric shape of the vessel itself. For 
instance, Tellep and Harper (Reference 4) employed a cylindrical con- 
tainer for analyzing the stratification of a liquid (water) receiving 
a constant heat flux from the surroundings. Ruder (Reference 5) pres- 
surized a contained fluid and discovered that the temperature profile 
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in the stratified layer had a shape similar to a Gaussian probability 
distribution function. 

Eckert and Jackson (Reference 6) do not deal specifically with 
stratification, but present an analytical study of a turbulent free 
convection boundary layer on a flat plate. The equations derived in 
their work provided the boundary layer temperature and velocity profile 
expressions used in the present study. 

Schwind and Vleit (Reference 7) present a very descriptive dis- 
cussion of the entire stratification process including the phenomenon 
of the reverse shear layer. Briefly, the reverse shear layer forms 
in a region between the boundary layer and the stratified section of 
the fluid core. This phenomenon occurs due to very complicated fluid 
particle motions. 

Unlike the side wall heating arrangement in the present investiga- 
tion, Duke (Reference 8) heated a vessel at the bottom. This vessel 
was cylindrical in shape and had a conical bottom. Despite the signifi- 
cant difference in design between this vessel and the one used for the 
present study, the temperature profiles in the respective stratified 
layers were found to be similar in shape. 

Bailey, et al. (Reference 9) studied the phenomenon of stratifica- 
tion of a cryogenic propellant. The analytical work was based on the 
approximate, integral method which was also used in the present study. 

The container was cylindrical in shape. The equations were essentially 
identical to those used by Tellep and Harper (Reference 4) . In this 
study, however, the analytical results were compared with actual data 
obtained from the propellant tanks of the Titan and Vanguard missiles. 
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Churchill (Reference 10) was not concerned with stratification, 
rather the significance here lies in the mathematical method used in 
the present study. The primary discussion was concerned with the 
attempts which have been made to predict natural convection phenomena. 
The study compares the classical methods of analysis, such as the 
approximate integral method employed in the present study, with the 
various new numerical techniques. 

A great deal of effort has been devoted to the phenomenon of 
stratification. Many of these previous studies were conducted with 
cryogenic fluids. These reports are listed under Additional References 
to provide a broader background for the present study. 
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DIMENSIONAL ANALYSIS 


In previous sections it has been explained that free convection 
heat transfer is the primary mechanism for energy transport in the 
stratification process. Free convection heat transfer has been thor- 
oughly studied by many investigators. The results of these investiga- 
tions disclose the parameters involved in the free convection heat 
transfer process, and consequently, the parameters imposed on the 
stratification process. In view of the knowledge of the free convection 
parameters, the stratified layer thickness (A) can be appropriately 
represented by 


A = f 


P, c, k, t, g, 3, y, 


<VV ’ H ’ 



( 1 ) 


The system is illustrated in Figure 3. 

The stratified layer thickness (A) is dependent on: 

(1) the elapsed time (t) 

(2) the physical properties of the fluid such as density (p) ; 
viscosity (y) ; conductivity (k) ; specific heat (c) ; and 
thfe thermal coefficient of expansion (3) 

(3) the gravitational field (g) 

(4) the heat transfer rate which is directly related to the 
temperature difference (T "T^) 

(5) the dimensions of the vessel, such as the height (H) ; 
the heated perimeter (P) ; and the area of fluid cross- 
section (A g ) . 

Equation (1) above expresses the stratified layer thickness (A) as a 



Adiabatic 


FIGURE 3. STRATIFICATION PARAMETERS 
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function of these parameters. Table 1 below lists these parameters 


with their corresponding fundamental dimensions of mass (M) ; 
length (L) ; and temperature (T) . 

Table 1. Fundamental Dimensions 

time (t) ; 

Symbol 

Parameter 

Dimensions 

A 

stratified layer thickness 

L 

k 

thermal conductivity 

-3 -1 
MLt T 1 

p 

density 

ML" 3 

c 

specific heat 

L 2 t’ 2 T" 1 

t 

time 

t 

g 

gravitational field 

Lt" 2 

6 

coefficient of expansion 

T' 1 

<vv 

wall to bulk liquid temperature difference 

T 

H 

height of liquid in the vessel 

L 

y 

dynamic viscosity 

ML _1 t 1 

P 

width of heated wall 

L 

A s 

cross-sectional area of the liquid 

L 2 


The dimensions associated with equation (1) can be represented by 


L = 



|ml 

rj 


3] b f L 


L 2 t' 2 T" J 


x ] c H d £ 
[ ML ' lt " 1 ] 1 H j H 


Lt 

k 


( 2 ) 


Since L alone exists on the left-hand side of equation (2), the sum of 
the exponents of the three remaining dimensions on the right-hand side 
must be zero. Accordingly, the sum of the exponents of L on the right- 
hand side of the equation must be one (1) . Listed below in Table 2 are 
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the four (4) dimensions and their corresponding exponents summed to the 
appropriate values as indicated above. The ensuing operations follow 
the method of Buckingham's n -Theorem. 


Dimension 

M 

L 

t 

T 


Table 2, Exponent Equations 
Exponents 

a + b + i 

a -3b +2c + e + h - i + j +2k 

-3a -2c + d -2e - i 

-a - c - f + g 


0 

1 

0 

0 


The four (4) equations contain eleven (11) unknowns; thus, four (4) of 
the unknown exponents can be solved in terms of the other seven (7). 

The exponents chosen are a, c, g, and h, respectively. These exponents 
are found in terms of the remaining exponents using the equations as 
given in Table 2. Equation (1) can now be written as 


A = (k) ■ b ' i tp) b (c) 3/Zb+i - e+l/2d (t) d (g) e (6) f 

(T -T, ) f>1/2b " e+1/2d (H) l +b ~ d+e 'j" 2 k (u )i f P )j (A ) k (3) 

v w b J s 


Rearrangement of this equation leads to 

)b 


f- [k^pc 372 (T w -T b ) 1/2 H 


1/2 _ T ) 1 / 2 h -1 
v w b J 


c -1 g (T -T, ) _1 H 
s ^ w b 


6 CT„-T b ) 


c v k 


PH 


1 

j 

( -2) 
A H z 

J 


s j 


( 4 ) 
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Equation (4) indicates that 



PC 372 (T w -T b ) 1/2 H 

j- 


: l72 t 


(T -TJ 
v w tr 

IT 


1/2 


£ 


(T ) 
v w b' 



(5) 


Some of the seven (7) n- factors which appear on the right-hand side 
of equation (5) can' be rearranged to give more familiar dimensionless 
groups. By carrying out the necessary manipulations, the expression for 
can be rewritten as 


FT 


= G 


[m 


pc 

» T 


» Z g8 CT W -T b ) H 3 


P> 




^W’ 


pc 3/2 (VV 1/2h P A s 
E • FT • J? 


( 6 ) 


The dimensionless groups in expression (6) are given descriptive names 
below. 


( 1 ) 

( 2 ) 

(3) 

(4) 

(5) 

( 6 ) 


jp dimensionless stratified layer 

D 

pp geometry of container 


-y, geometry of container 
FT 

, dimensionless time 


vH7g 


gp 2 e (T w -T b ) H 3 


, Grashof Number 


u 

Prandtl Number 


thickness 
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( 7 ) 

( 8 ) 


rT ^ - y , ratio of potential energy to thermal energy 

cl W 


pc 3/2 (T w -T b ) 1/2 H 


, dimensionless diffusion of thermal energy 


The above results show 'that no less than eight (8) dimensionless 
groups are pertinent to the Stratification process. This information 
serves as a guide for the analysis to follow. 



ANALYSIS 


There are two major objectives of the analysis presented in this 
section. The first is to obtain a relation between the stratified layer 
thickness (A) and elapsed time (t) . The second is to determine the 
temperature distribution in thte stratified layer and, simultaneously, 
to obtain the liquid surface temperature (T s ) variation with time. The 
fluid is considered to have ari' initially uniform temperature (T^) . The 
various parameters used in the analysis are shown in Figure 4. 

The first portion of the ‘analysis involves the derivation of an 
expression for the growth of the stratified layer based on the mass 
flow rate (m^) in the boundary layers. From this foundation and from 
a consideration of the internal energy of the stratified layer, it is 
possible to determine the approximate temperature distribution and the 
rate of increase of the liquid surface temperature (T ) . 

Certain assumptions are made for the purpose of promoting a straight 
forward step-by-step analysis resulting in the determination of these 
expressions. This procedure leads to analytical predictions which are 
compared with experimental measurements . The assumptions made in the 
analysis are the following: 

1. The heat transfer rate from the heated walls to the contained 
■ fluid remains constant throughout the transient stratification 
period, and equal to the actual rate existing just after the 
process is initiated. In reality, the heat transfer rate 
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Legend: 

q - Heat Transfer through Walls 
6(x) - Boundary Layer Thickness 
A - Stratified Layer Thickness 
T w - Inside Wall Temperature 
T g - Surface Temperature 
T^ - Bulk Temperature 

- Mass Flow Rate in Boundary Layers 



FIGURE 4. DESCRIPTION OF PARAMETERS USED IN ANALYSIS 
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decreases with time to some extent, since the heat flux from 
the portion of the heated walls adjacent to the stratified 
layer is attenuated by the increased fluid temperature levels 
in the upper portion of the tank. Consequently, the rate of 
increase of the surface temperature (T ) as determined by this 
analysis is somewhat higher than that which actually occurs. 

2. There is no heat transfer from the boundary layers or from 
the stratified layer to the bulk fluid below the stratified 
layer; therefore, the bulk temperature (T, ) remains constant 

.1 

throughout the stratification process. 

3. The convection heat transfer coefficient (h) is calculated 
according to (Reference 11) : 

h = 0.021 (GrPr) 2/5 (1) 

for a turbulent boundary layer (GrPr > 10 ^). However, when 

8 10 

the GrPr product is in the neighborhood of 10 - 10 , the 

heat transfer coefficient is determined by using the following 
correlation (Referende 11) : 

k 

h = 0.12 ~ (GrPr) 1/3 (2) 

The corresponding expression for a laminar boundary layer is 
not presented here, since it is assumed that the boundary 
layers existing in typical stratification processes are either 
in the transitional or fully turbulent range. (This contention 
is illustrated in a subsequent section of this report which 
deals with the experimental results.) 
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4. All fluid properties are assumed to be constant. 

5. The temperature profile in the stratified layer is assumed to 
be well approximated by a parabola (see Figure 4) . 

6. The boundary layers are very thin such that the boundary layer 
response can be assumed as instantaneous compared with the rest 

of the system. Thus, it is assumed that steady state velocity and 
temperature boundary layer profiles apply. Hence, the thermal 
inertia of the boundary layers is neglected. 

The statements 1 through 6 above complete the list of assumptions and 


conditions incorporated into the analysis. 

The thermal conditions discussed under assumption 2 require some 


further explanation. The heated region at the top of the container 
in reality extends slightly below the stratified layer as a result 
of heat conduction between the stratified layer and the colder 
bulk fluid below it. However, in order to simplify the analytical 
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description of the stratification process, it is assumed that the 
heat transfer occurring between the stratified layer and the bulk 
* fluid is negligible. This implies that the thickness of the heated 
layer at the top is essentially equal to that of the stratified layer . 
This is a reasonable assumption since transient heat conduction effects 
can be shown to be small for time periods which are typical of the 
stratification process, as mentioned previously. 

The goal of the present analysis is to formulate a simple mathe- 
matical model of the stratification phenomenon by considering only the 
dominant physical processes. Thus, the resulting analytical predictions 
are necessarily approximate and must be regarded as qualitative, pri- 
marily due to the fact that the extension of the heated region below 
the stratified layer is not treated. Nevertheless, the analysis does 
yield approximate solutions which provide the proper trends and orders 
of magnitude for the important parameters of interest. This fact was 
established on the basis of experimental measurements to be presented 
later in this report. These measurements also indicate that the growth 
of this heated region does not appreciably influence the rate of growth 
of the stratified layer until far into the transient period when the 
heated portion of the bulk liquid extends nearly to the bottom of the 
container. 

Refer now to Figure 4 and consider the boundary layer to be turbu- 
lent. Eckert and Jackson (Reference 6) give approximations to the 
velocity (u) and temperature (6) profiles for a turbulent boundary 
layer as 


u = u l (y/6) 1/7 (l-y/6) 4 


w 


1 - (y/6) 


1/71 


(3) 


0 = 6 


(4) 
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where 

= velocity outside of the boundary layer of the comparable 
forced convection flow 

e = deviation of temperature within boundary layer from the bulk 
fluid temperature = (T-T b ) 

0 ^ = temperature difference between the wall and the bulk fluid = 

<VV 

u = velocity of the fluid within the free convection boundary 
layer 

The expression for u^ is also obtained from Reference 6 as 

_ 1/2 

Ul (x) = 1.185 (v/x) Gr 1/2 jl + 0.494 (Pr) 2/3 | (5) 

where 

x = coordinate measured upward from the container bottom 
6 (x) = boundary layer thickness 
Also from Reference 6 , 6 (x) is given by 

6(x) = 0.565x'(Gr)" 1/10 (Pr)' 8/15 1 + 0.494 (Pr) 2/3 J Vl (6) 

The mass flow rate in the two boundary layers, at a distance x 
from the vessel bottom, is obtained by integrating the velocity (u) 
over the boundary layer thickness. Thus, 

r6(x) 

m K1 (x) = 2P pudy (7) 

01 Jo 

where P is the width of each heated wall. Since the value of the fluid 
density (p) is taken to be constant, the expression for the mass flow 


rate becomes 
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f<5 (x) 

m bl (x) = 2Pp udy (8) 

The relation for the velocity as a function of distance from the wall, 

equation (3) , is substituted for u in equation (8) . The result is , 

■6 (x) 1/7 A 

I^ 1 (x) = 2p p J o u x (x) [y/fiCx)] 177 [1 - y/<S (x) ] dy (9) 

The variable of integration (y) is now changed 'to y/6(x) by rewriting 
equation (9) as 

m bl (x) = 2Pp6 (x) u-^x) [y/6(x)] 1/y [1 - y/6(x)] 4 d[y/6(x)] (10) 

Integration of this expression yields 

m bl (x) = 2Pp [0.145 6(x) u^x)] (11) 

The previously given expressions for u-^(x) and 6(x), equations (5) and 
(6), are now substituted into the above equation, and the final expres- 
sion for the mass flow rate is obtained as 

m bl (x) = 2Pp< 0.097 v (Gr) 2/5 (Pr) 8/15 [l + .494 (Pr) 2/3 ' j (12) 

Since the growth of the stratified layer comprises the first objective 
of this analysis, the mass flow rate (m^) is evaluated at its location 
of entrance to the stratified layer. Specifically, the mass flow rate 
(m bl ) is calculated at the lower boundary of the stratified layer where 
x = H - A. Thus, equation (12) becomes 
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lV x) ] x . h . ; 2Pp °- 


097v 


gg (T w -T b ) (H-A)- 


2/5 

(PrJ 


8/15 


■2/S 


Q + . 494 (Pr) 2/ ^j | (13) 

which gives the mass flow rate across the lower boundary of the strati- 
fied layer- Notice that equation (13) can be written in simplified 
form as 

m K1 (x)l = Constant (H-A)^^ •= C-, (H-A)^ (14) 

- J x » H - A 


where C^ = 0.194 Pp v 


**<V T b ) 


2/5 


(Pr) 8/1S 


1 + „494(Pr) 


) 2/3 ]- 


2/5 

lbm 


hr-ft 


Vs 


Another expression containing the boundary layer mass flow rate (m^) 
can be developed by applying the law of conservation of mass to the 
stratified layer. This expression is 


V00 


1 at (pA s A} = pA s 3? 


(15) 


-> x = H - A 
Combining equations (14) and (15) results in 

» dA t rj , 1 6 / 5 

i>A s at c i lH 01 

Equation (16) can be integrated to determine the relation for A as a 
function of time (t). Thus, separating variables leads to 


(16) 


dA 


C^dt 


(H-A ) 6/5 pA 


(17) 


5(H-a) 1/5 - — ♦ K. 

PA, 


( 18 ) 


Integration produces 
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where K 1 is a constant of integration. is determined by applying 
the boundary condition that A = 0 at t = 0. From equation (18) the 
result is 


K x - 5H~ 1/5 (19) 

Consequently, equation (18) can be rewritten as 

5CH-A)" 1 / 5 = — + 5H' 1 / 5 (20) 

pA„ 


The stratified layer thickness (A) can be obtained as a function of 
time (t) through the rearrangement of this expression. Thus, 



1 



SpA s 


( 21 ) 


The constant factor can now be replaced by its corresponding para* 
metric expression as obtained from equation (13) . This operation yields 
i 


A 

H 


1- 

1 + .0388 v 

r P ' 

a 

(Pr) 8/lS 

( T w' T b^ 

1 


A 

t 


2 

v _ 


2/5 




1 + .494 


(Pr) 2/3 


-2/5 


-5 




( 22 ) 


Equations (21) and (22) express the growth of the stratified layer 
thickness (A) as a function of time (t). It remains at this point in 
the analysis to determine the variation of the surface temperature (T ') 
as a function of time (t)« The surface temperature (T ) can be inserted 
into the expression which equates the rate of convective heat transfer from 
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the heated walls to the change in internal energy occurring within the 
stratified layer „ This expression is 


1 - ’V J 0 < T - T b’ dz 


where 


q = heat transfer rate in Btu/hr 

2 

A = cross-sectional area of fluid in ft 
s 

t = elapsed time in sec 
c = specific heat in ^^ U ep - 

A = stratified layer thickness in ft 

z = coordinate measured from the liquid surface downward 
Equation (23) is now rewritten as 









at \ 

ii 

i ( T "V A 1 
s b J 0 

- T , T n 


f ' 

z 


s-L 



(23) 


(24) 


\ / 

In order to carry out the integration of equation (24) , the quan- 


tity 


T-T, 


T T, 
s b 


must be known as a function of z/A. At this point assumption 5 is 
employed. Accordingly, the approximation for the dimensionless tempera- 
ture profile in the stratified layer is written as 

T-T 

Y~Y- “ U • z/a ] 2 (25) 

s b 


This approximation has been used previously by Tellep and Harper 
(Reference 4), who offered considerable evidence of its validity for 



27 


their investigation. The justification for using this procedure in 
the present development is based on the experimental temperature mea- 
surements by Tellep and Harper as well as those of the present study. 
Tellep and Harper found that use of this expression for the temperature 
profile gave rise to a predicted rate of growth for the stratified layer 
which compared very well with their experimental results. However, 
there was a sizeable deviation between the analytical prediction and 
the measured data when the stratified layer approached the container 
bottom. Tellep and Harper's results are compared with those of this 
study in a following section of this report. The insertion of equation 
(25) into equation (24), followed by evaluation of the integral, pro- 
duces 

q = 1/3 pA s c ^ [(T s -T b ) A] (26) 

An exact integration of equation (26) is not achievable at this 
stage of the analysis, since the time dependency of the left-hand side 
(q) is not known. Here it is assumed that the heat transfer through 
the heated walls remains constant throughout the entire transient period, 
and equal to the actual rate occurring immediately after initiation of 
the process (assumption 1) . Thus , 

q = hA (T -T, ) = Constant (27) 

n w b 

where 

q = heat transfer rate in Btu/hr 

h = average convection heat transfer coefficient at the start 

of the process m -r— - 

hr-ft °F 
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A = total wall area from which heat is transferred to the 
contained fluid, 2PH in ft^ 

P = width of the wall in ft 
H = total depth of the fluid in ft 
As mentioned previously in the discussion of assumption 1, the 
actual heat transfer rate can be expected to decrease with time to some- 
what lower values as stratification proceeds. Therefore, the predicted 
variation of surface temperature (T ) with time (t) can be expected to 
be above that which occurs in reality. Insertion of equation (27) into 
equation (26) yields 

q = fiA (T w -T b ) = Constant = 1/3 pA g c ^ [ (T s _ T b ) A] ( 28 ) 


or 


^ [(T s -T b ) A] = ^ = K (constant) (29) 

Integration of equation (29) yields 


(T -TJ A = Kt + E (30) 

s b 

where E is a constant of integration. Consideration of the two initial 
conditions, A = 0 and (T -T^) = 0 at t = 0, leads to the conclusion 
that the constant E must equal zero also. Equation (30) is now solved 
for the temperature difference with A replaced by its equivalent expres- 
sion given in equation (21). Thus, 


CTe-TJ = 


Kt 


Kt 




C 1 H 1 / 5 t 

-5' 


1 - 


> 




. 


(31) 
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Substituting the expression for K from equation (29) yields 


CVV 


(T w -T b ) t 



' 

C.H 1 ^ 

-5 

0 A s Hc < 

1 - 

1 + 5pA_ 
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(32) 


The combination of parameters, pA s H, in the denominator of equation (32) 
is actually the total mass of liquid in the container. Equation (32) 
can be nondimens ionali zed by dividing through by the temperature differ- 
ence (T*T b ) which yields 


T -T, 
s b 


6hPt 


w 


»A s C 


a 

C,H 1//5 t 

1 + ^ . 

■ 5 1 



1 

— — 



(33) 


where the heat transfer area, A, has been replaced by 2PH, A final 
expression for the temperature difference is obtained by substituting 
the full expression for the constant from equation (13). Thus, 


T 

s 


T 

w 




( 34 ) 


where A , the liquid cross-sectional area, has been replaced by the 
product of P and L which are the heated and unheated wall widths, 
respectively (see Figure 3).. 
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The expressions developed above for the predicted growth of the 
stratified layer and the increase in the surface temperature during 
the stratification process, equations (22) and (34) respectively, are 
plotted and compared with the results of Tellep and Harper (Reference 4), 
as well as experimental measurements from this study, in Figures 10 
through 13. The numerical calculations used for these comparisons are 
presented in the Appendix of this report. 
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EXPERIMENTAL APPARATUS AND PROCEDURE 


The vessel used in the experimental program (Figure 5) was made 
primarily of plexiglass. Two opposite sides each had a one-eighth (1/8) 
inch thick aluminum sheet inserted approximately three-quarters (3/4) 
of an inch inward from the plexiglass walls. The aluminum sheet together 
with the adjacent plexiglass wall created a channel which extended over 
the entire height and width of the vessel. These channels provided 
a path through which hot water flowed and consequently heated the alu- 
minum walls. The top of the vessel was constructed from a combination 
of wood and Styrofoam which had a high thermal resistance, thereby 
minimizing the loss of heat and evaporation at the liquid surface. 
Plexiglass ribs were attached to the channel side of the aluminum sheet. 
The function of the ribs was to maintain a uniform flow in the vertical 
direction throughout the channel. Water was used as the test liquid 
in all of the experiments performed. 

For the purpose of fluid temperature data acquisition, a thin 
wooden post was positioned vertically in the center of the vessel as 
depicted in Figure 5, The post was attached to the bottom with an 
adhesive. Chromel- alumel thermocouples were attached at specific inter- 
vals along the height of the post. With the interior of the vessel 
filled to a depth of twelve (12) inches with water, the vessel width 
was six (6) inches and the heated wall width was twelve (12) inches. 

The uppermost thermocouple was located precisely at the water surface. 
The remaining thermocouples were located at the following depths below 


the surface: 
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FIGURE 5. VESSEL WITH CENTERLINE THERMOCOUPLE POST AND FLOW SEPARATORS 
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one-half (0.5) inch, one (1.0) inch, two (2.0) inches, six (6.0) inches, 
and eleven (11.0) inches. The reference junction was located in an 
ice bath which was maintained at 32 °F. A digital voltmeter equipped 
with a tape printout system completed the thermocouple circuit. This 
instrument is shown in Figure 6. The digital voltmeter (DVM) had the 
capability of reading ten (10) inputs sequentially and cyclically . 

Each printout consisted of a voltage reading plus the number signifying 
the specific thermocouple being read. A complete reading and printing 
operation for each thermocouple was performed within an interval of 
one second. Other thermocouples were placed at the top and bottom of 
the inside and outside of the aluminum walls. These temperatures were 
also recorded by the DVM and provided a record for the wall temperatures 
throughout each experiment. 

The flow through the wall channels originated from a tank of hot 
water which was located above the vessel. Flow separators, shown in 
Figure 5, had to be inserted in the hot water flow stream between the 
source and the vessel channels to insure a proper uniform flow through 
the channels. When a valve was suddenly opened to allow the hot water 
to flow through the two side channels, the thin aluminum sheets were 
heated very rapidly to a temperature just slightly below that of the 
hot water. This provided a good approximation to a step change in 
wall temperature as far as the interior fluid was concerned. 

Experiments were conducted in order to obtain measurements of 
some of the important variables in this study, namely, the liquid sur- 
face temperature (T ) , the transient centerline temperature distribu- 
tion, and the stratified layer thickness (A) . 



FIGURE 6. EXPERIMENTAL APPARATUS COMPLETE WITH AUXILIARY EQUIPMENT 
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FIGURE 6 




Each experiment was initiated by a valve which allowed hot water 
to flow through the separators and through the vessel channels along 
the two aluminum walls. Heat transfer through the aluminum walls cre- 
ated free convection boundary layers along the interior walls, thus 
initiating the stratification process. The hot water source was thermo- 
statically controlled and the hot water temperature was measured at 
the entrance to each channel. 

The centerline post, shown in Figure 5, held six f6) thermo- 
couples, including the one at the surface. This surface thermocouple 
was used to obtain measurements of T . The DVM system printed out 
values of the surface temperature, as well as the other temperatures 
along the vertical centerline, as indicated by the various thermocouples. 
The data was converted from millivolts to degrees Fahrenheit. Each 
measurement was made dimensionless by subtracting the bulk temperature 
(T,) from it, and then dividing this result by the temperature differ- 
ence (T -T, ), Plots of the resulting measurements, obtained with 
several of the thermocouples, are given in the next section of this 
report. The increase in the stratified layer thickness (A) was deter- 
mined by the initial change in each thermocouple reading. Thus, the 
downward propagation of the stratified layer, as observed experimentally, 
is presented.. 

The stratification process was also observed by means of a blue 
dye. The dye mixture consisted of iodine, water, and starch and had 
a specific gravity of 1.015. The dye was introduced into the bulk 
fluid, where time was allowed for it to settle to the bottom and to 
reach the same temperature as that of the bulk fluid (T^j , before each 
observation was initiated. Note that the dye was slightly heavier 
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than water. This condition was chosen so that the dye would descend 
to the bottom prior to each visual test, so that it would clearly show 
the upward movement of the fluid once the transient heating process 
was initiated. The dye settled on the bottom with a thickness somewhat 
less than an inch. Consequently, the dye contacted the boundary layer 
near its point of initiation on the wall. Unfortunately, at this loca- 
tion the boundary layer was laminar with a very small velocity and was 
unable to raise the dye for the lower temperature differences. Only 
for the largest temperature differences used (T -T^) was there a suffi- 
cient flow which was capable of raising the dye. This condition implies 
that the dye flow did not coincide with the front of the stratification 
process. In fact, in comparing the centerline thermocouple readings 
with the dye flow, it became clear that the dye followed significantly 
behind the stratification front. Nevertheless, the dye flow provided 
a clear means by which the formation of the stratified layer, as illus- 
trated in Figure 2, could be observed. 

The dye flow was observed with a television camera and recorded 
on a video tape recorder. A clock (Figure 6) was also in the video 
picture. It provided an account of the elapsed time as the dye flowed 
up the boundary layer and into the stratified layer. Through playback 
of the video tape, the mixing occurring within the stratified layer 
could be observed as many times as desired. Later in the experimental 
program a motion picture camera was substituted for the television 
system. This was done to provide a lasting account of the process. 

Plus X negative film was used with the camera which was set at sixteen 
(16) frames per second. 
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The two hot water temperatures used were 140 and 180 °F. Thermo- 
couples were positioned at the bottom and top on both the inside and 
outside of the aluminum walls. With these thermocouples the actual 
wall temperature gradients were recorded throughout the experiment. 

The minimum elapsed time for all of the experiments was eight (8) 
minutes. This allowed time for each of the six (6) thermocouples to 
record several temperature values which provided enough plotted data 
points so that fairly smooth curves could be drawn through them. 



RESULTS 


Typical data obtained with the thermocouples , employed as described 
in the preceding section, are presented in Figures 7 and 8 for the 
centerline temperatures in the fluid. The surface temperature (T ) 
lags momentarily, rises rapidly, and then increases at a smoothly 
decreasing rate while asymptotically approaching the Value of the inside 
wall temperature (T ) . Experiments were conducted at two different 
values of the temperature difference, arid its effect can be evaluated 
by comparing Figures 7 and 8. A comparison shows that the higher tem- 
perature difference causes the centerline temperature's to increase 
more rapidly to an asymptotic value. The higher temperature difference 
(Figure 8) shows a greater increase in the centerline temperatures 

rtf 

at the depths of one (1.0) and six (6.0) inches. This latter result 
indicates the increase in' the rate of growth of the stratified layer 
for the higher temperature difference compared to the growth existing 
at the lower temperature difference. Figure 9 contains plots of the 
centerline temperature as a function of depth at three different times. 
The total depth of the fluid was twelve (12) inches. The data show 
that it took approximately two (2) minutes before the centerline tem- 
perature began to increase at three-quarters (3/4) of the total depth. 

The analytical predictions formulated in this study for the growth 
of the stratified layer thickness (equation 22) and the rate of increase 
of the surface temperature (equation 34) are compared with typical 
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experimental data and with Tellep and Harper's analytical results in 
Figures 10 through 13. The numerical calculations associated with the 
analytical expressions involved are contained in the Appendix. The 
analytical prediction for the growth of the stratified layer thickness, 
as displayed in Figures 10 and 11, compares quite favorably with the 
experimental measurements. The data points shown were obtained by 
noting the time of the initial increase in temperature for each of 
the centerline thermocouples which were located at specific depths 
below the surface. The analytical prediction lies below the data during 
the first thirty (30) seconds of the stratification process and is 
somewhat higher thereafter. The analytical prediction of Tellep and 
Harper (Reference 4) falls below both the experimental data and the 
analytical prediction of the present study over the entire range of 
elapsed time shown in Figures 10 and 11. However, the general shape 
of the curve, including the variation of its slope with time, compares 
well with both the experimental data and the analytical prediction 
obtained in this study. The measured results show that about 801 of 
the twelve (12) inch depth of fluid was stratified after only two (2) 
minutes had elapsed. The curves of Figure 11 relate the stratified layer 
thickness to dimensionless time, and therefore, to fluid and container 
parameters . 

An analytical prediction of the surface temperature (T ) increase 
during the stratification process is shown for comparison with some 
typical measurements in Figures 12 and 13. The order of magnitude of 
the prediction correlates very well with the experimental data during 
the time in which 80% of the contained fluid is being stratified. 

However, it can be seen that the predicted surface temperature variation 
has a discontinuity at t=0. Obviously this is not realistic, as indicated 
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by the smoothly increasing trend observed experimentally , and is due to 
the neglect of the boundary layer inertia (assumption 6) in the analysis. 

It can also be seen that after two (2) minutes have elapsed, the analytical 
prediction rapidly diverges from the experimental curve. This result can be 
attributed largely to the assumptions used in the derivation of equation 
(34) , particularly assumption 1 which implies a constant heat transfer 
rate. This becomes unrealistic as the stratified layer approaches the 
bottom of the container. Note that equation (34) states the surface 
temperature rate of increase is directly affected by the growth of the 
stratified layer thickness during the entire stratification process. 

The experimental data, on the other hand, show that after about 80% 
stratification has occurred, the surface temperature changes very little 
with further elapsed time. 

Tellep and Harper's analytical prediction for the surface temperature 
increase is also shown in Figure 12 and 13. It can be seen that their 
prediction does not correlate well with the experimental data from this 
work in the early portion of the transient region, but does become 
approximately coincident with the measurements in the 80 to 90% stratification 
range. Beyond this range their prediction lies well above the experimental 
data points, but the divergence from the experimental results in the latter 
portion of the transient period is only about one half as large as the 
deviation exhibited by the analytical prediction of the present study. 

It should be noted that the analytical predictions of the surface temperatures 
should logically lie above the measured data for large times, since both 
analyses incorporate the assumption of a constant heat transfer rate, 
whereas, the measured data were obtained with essentially a constant wall 
temperature condition leading to an actual heat transfer rate which 
eventually decreases considerably as the elapsed time increases. 
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FIGURE 8. OBSERVED TEMPERATURE CHANGE AT DIFFERENT DEPTHS WITHIN THE STRATIFIED LAYER (T w -T b = 56.2°F) 
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FIGURE 9. OBSERVED TEMPERATURE PROFILES ALONG THE VERTICAL CENTERLINE (T -T 
















CONCLUSIONS 


This study has involved an analytical and experimental investiga- 
tion of the transient stratification of water contained in a rectangular 
vessel. It is concluded that the stratification of a contained fluid 
with side wall heating due to a constant temperature can be approxi- 
mately predicted by means of an analysis based on a simplified physical 
model. A summary of the pertinent results is given below: 

1. An analysis based on fundamental heat transfer equations 
produced an expression which predicted approximately the 
growth of the stratified layer. This was confirmed experi- 
mentally by observing the times at which the various center- 
line temperatures first began to increase. 

2. The same analysis yielded an expression for the increase 
of the surface temperature which compared favorably with 
experimental data during the early portion of the transient 
period, but which was found to overpredict the surface tem- 
perature during the latter portion of the transient period. 

3. The dimensionless results indicate the manner in which the 
fluid properties and vessel geometry affect the stratification 


process. 
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APPENDIX 


NUMERICAL CALCULATIONS 


In this section the analytical expressions will be utilized for 
a specific condition. The temperature difference (T^-T^) to be inves- 
tigated is 

(T -T,) = 135.9 °F - 79.7 °F = 56.2 °F. 

v W b' r 

The parameters will be evaluated at the average temperature 

T +T, 

f = = 107 °F . 

The tables in Reference 11 provide the following values: 

Pr = 4.52, 6 = 0.219 x 10" 3 /°F, v = 0.692 x 10 -5 ft 2 /sec. 


The known parameters are, 

P = 1 ft, g = 32.2 ft/sec 2 , A = 1/2 ft 2 , L = 1/2 ft, H = 1 ft. 


Let’s recall equation (22), 


A 

FT 


( 1 + .0338 v 

' P' 

m sni [ 

T 


s . 

- 


l 2/5 


[l + .494 (Pr) 2/3 J / H 1/5 t 
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which after termwise exponentiation becomes. 


A. i . 

H L 


1 + .0388 


2 

7t 


0.692 x 10 


-5 ft2 
sec 


(2.24) L 926 * 10 4 /ft 6/5 J 
-5 


[1 + .494 (2.738)]' 2/5 (ft 1/5 ) t(hr) } 
This expression reduces to 


H = 1 " \ 1 + 28,44 Hr / 

^ values can be obtained by substituting time in hours for t. The 
following values result. 


t (hr) 

28.44 t 

1 + 28.44 t 

(1 + 28.44 t)' 5 

A/H 

0 

0 

1 

1 

0 

.001 

.028 

1.028 

.872 

0.128 

.005 

.142 

1.142 

.516 

0.484 

.010 

.284 

1.284 

.288 

0.712 

.015 

.427 

1.427 

.170 

0.830 

.020 

.569 

1.569 

.106 

0.894 

.040 

1.138 

2.138 

.023 

0.977 

.050 

1.422 

2.422 

.012 

0.988 

.060 

1.706 

2.706 

.007 

0.993 

.080 

2.275 

3.275 

.003 

0.997 

.100 

2.844 

3.844 

.001 

0.999 


The values for ^ from this table are plotted and compared with Tellep 
and Harper’s results, and experimental data from this work, in Figure 


10. 
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The calculations continue by employing equation (34) to obtain the 
analytical rate of increase of the surface temperature (T s ) for the same 
specific condition used to obtain (A) in the preceding calculations. 

The surface temperature increase is predicted by equation (34) , 


V T b 

T -T, 
w b 


6h 


IpEc 




1 + .0388 v 


P' 


(Pr) 


8/15 


^vv 
2 


2/5 


1 + .494(Pr) 


2/3 


-2/5 


H^t 


-5 


The heat transfer coefficient h will be calculated using the expression 
for turbulent flow given in the analysis as equation (1). Thus, 


h = 0.021 (GrPr) 2/5 


where 


Pr = 4.52 


gB(T -T b )H° 

Gr = = 825 x 10 ' 


k£ = 0.363 


Btu 


hr-ft F 


and 


h = .021 (.363) 
Btu 


(825 x HT) (4.52) 


-l 2/5 


Btu 


hr-ft 2 °F 


fi = 129.0 


hr-ft 2 °F 
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The constant term in the numerator can be calculated, 


= 25.048/ hr 


which is substituted into equation (34) along with A/H for the denomi 
nator c Thus, 


T -T, 
s b 

w b 


25.048 t 
hr 
A 
H 


The values for T can be found by substituting time in hours for (t) 
and the corresponding values for ^ from the preceding table. 


(hr) 

A/H 

T -T, 
s b 

T °F 
s 

IR x 10 ‘ 

.001 

.128 

10.96 

90.66 

.118 

.005 

.484 

14.54 

94.24 

.593 

.010 

.712 

19.77 

99.47 

1.186 

.015 

.830 

25.44 

105.14 

1.779 

.020 

.894 

31.49 

111.19 

2.372 

.040 

.977 

57.62 

137.32 

4.744 

.050 

.988 

71.26 

1 150.96 

5.930 

.060 

.993 

85.04 

164.74 

7.116 

.080 

.997 

112.92 

192.62 

9.488 

.100 

.999 

140.94 

220.64 

11.860 


These data are plotted and compared in Figure 12. 

Since Tellep and Harper's work (Reference 4) played a very important 
part in the analysis, the results, both analytical and experimental, 
were compared with their data in Figures 10 through 13. It must be 
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understood that their work was conducted in a cylindrical container as 
compared to the rectangular container used in this work. This conflict 
was considered and it was reasoned that the radius of the circumscribed 
circle of a square using P = 1 ft as a side would be an approximation 
to the radius of the cylindrical container (R) . Therefore, since the 
length of a radius of a circumscribing circle is 3/4 of the length of 
a side, the following relation was used: R = 3/4, P = 1 ft. This 
value of R was substituted into Tellep and Harper's turbulent prediction, 

namely, 




= 1 - 


1 + .082 (H/R) (Gr*/Pr ) 2/7 <t> 


i-7 


where 


<P T 

n 


geqH 4 

Gr* = (q = wall heat flux) 

k£v 

The wall heat flux (q) was determined using the temperature difference 
(T w "Tb) incor P orate(i in the previous calculations. Thus, 

q = h (T -T, ) = constant (Reference 4) 
n w b 

Btu 


and 


q = (129) (56.2 U F) = 7249.8 

Gr* = 293.8 x 10 10 , Pr = 4.52 
<}> = .025 (t hr) 


hr- ft 


Values for jjj can now be found by the following procedure: 
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tlhli 

A/H 

0 

0 

.0032 

0.110 

.0064 

0.215 

.0129 

0.410 

.0193 

0.560 

.025 7 

0.670 

.0515 

0.870 

.0772 

0.940 

.1030 

0.975 


These values are plotted in Figure 10. 

A comparison was also made with Tellep and Harper's surface tempera- 
ture increase results. Their turbulent prediction is, 

T s -T b 2 (H/R) $ 

q H/k f I Pr (1 - [A/H] } 

where I is an energy integral representing the average temperature rise 
in the stratified layer. They suggest a value of 0.33 which was used 
to plot T . The following table gives the values used to plot T vs 
time in Figure 12. 
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t(hr) 

T -T, 
s b 

T °F 
s 

at . 

ra x 10 

0 

0 

79.7 

0 

.0064 

23.97 

103.7 

.76 

.0129 

29.96 

109.7 

1.53 

.0193 

31.95 

111.7 

2.29 

.0257 

33.95 

113.7 

3.05 

.0386 

43.94 

123.6 

4.58 

.0515 

51.92 

131.6 

6.11 

.0772 

71.90 

151.6 

9.16 

.1030 

93.86 

173.6 

12.22 
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AN INTERFEROMETRIC STUDY OF THERMAL STRATIFICATION 
BY NATURAL CONVECTION IN A CONTAINED LIQUID 1 


1 This section of Part V of the Final Report consists of a paper 
which was published in Warme- und Stof fubertragung , Vol. 4, 
September 25, 1968, pages 1-5. 



AN INTERFEROMETRIC STUDY OF THEFMAL 
STRATIFICATION BY NATURAL CONVECTION 
IN A CONTAINED LIQUID 

R. I. Haug, W. H. Stevenson, E. R. F. Winter 
Lafayette, Indiana, U.S.A. 

Abstract Transient laminar free convection to water in a rectangular 
container was investigated using a Mach-Zehnder interferometer . The 
experimentally observed temperature profiles are compared with analytically 

O 

predicted profiles based on a computer program developed by Barakat. Good 
agreement between analysis and experiment was obtained. In addition 
qualitative observations of the conditions in the container during sudden 
discharge of the liquid are reported. 

Zusantnenfassung Untersucht wurde der unstetige , laminare Wanteubergang 
durch freie Konvektion in Wasser in einem rechteckigen Behalter mittels 
eines Zehnder-Mach Interferometers. Die im Experiment beobachteten Tem- 
peraturprofile wurden mit den analytisch voraus ges agten Profilen verglichen, 
die mit einem Kamputerprogranin, das von Barakat entwickelt wurde, erzielt 
wurden. Gute Uebereinstimmung zwischen Analyse und Experiment wurde 
erreicht. Ausserdem werden qualitative Beobachtungen der Vorgange im 
Behalter wahrend des plotzlichen Leerens der strati fizierten Behalter- 
flussigkeit berichtet. 


INTRODUCTION 

Daring the past decade the process of thermal stratification by 
natural convection heat transfer to contained liquids has been the subject 
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of numerous analytical and experimental investigations. A review article 
by Clark [1] covered research done in this area through 1965. More recent 
analytical and experimental work is reported by Barakat [2] and Barakat 
and Clark [3], [4]. These latter studies are of particular interest here, 
since a complete numerical solution to the momentum and energy equations 
was obtained for the laminar flow case without invoking the boundary layer 
assumptions . The predicted temperature field near the heated vertical 
walls was found to differ markedly from that developed in an unconfined 
fluid. A prime objective of the present investigation was to experimentally 
verify this behavior using a Mach-Zehnder interferometer. This experi- 
mental technique also allowed certain other aspects of the stratification 
process to be observed qualitatively, including the temperature distri- 
bution in the liquid during discharge fran the container. 

Since the stratification phenomenon is so well documented, only a 
brief discussion of important considerations in the present study is 
necessary here. The governing equations for this problem are the x and y 
momentum equation, the continuity equation, and the energy equation. A 
computer solution of these equations was obtained by Barakat after first 
transforming the momentum and continuity equations into dimensionless 
expressions in terms of the stream function and the vorticity. Details 
of this transformation may be found in the publications by Barakat and 
Clark [2] , [3] , and [4] . The computer program given in [2] was modified 
slightly and used by the authors as a basis for comparison of experimental 
results with analytical predictions. 

To aid in the explanation of essential features of the problem the 
important assumptions made in the analysis are presented here. 
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Hie following thermal boundary conditions were used: 

(1) Hie bottom surface of the container and the free surface of 
the liquid are perfectly insulated. 

(2) Hie liquid in the container is initially at a uniform tem- 
perature T . 

(3) Starting at time zero the wall temperature T w (assumed to 
be uniform with location) is increased in some prescribed 
manner to a final steady value, T w (°°) . 

Hie boundary and initial conditions for the momentum equation are 
listed below: 

(1) No slip exists at the wall. 

(2) Hie shear at the free surface is zero. 

(3) Hie velocity is initially everywhere zero in the container. 

In addition to the above conditions , several other simplifying 

assumptions were made. Hie dissipation term in the energy equation was 
assumed negligible. Flow was assumed to be laminar and incompressible 
at all times and the process was taken to be two-dimensional. Also all 
fluid properties were assumed constant with the exception of the density 
in the body force term of the momentum equation which was allowed to vary 
with temperature . 

EXPERIMENTAL APPARATUS 

A Mach-Zehnder interferometer with a helium-neon laser source was 
employed to obtain data cn temperature profiles in the contained liquid. 
Hie use of a laser eliminated the need for a compensating chamber in the 
reference path of the interferometer. (A good discussion of the use of 
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a compensating chanter is presented by Goldstein and Eckert [5] , who used 
an interferometer to study transient free convection from a vertical plate 
in water. Another interferometric study pertinent to the present exper- 
iment was the investigation of laminar free convection heat transfer in a 
liquid with internal heat generation carried out by Novotny and Eckert [6] .) 

Figure 1 shows the test apparatus used in this study. It consisted 
of two independent units - a Plexiglass container with optical glass end 
windows to hold the liquid and a separate heating section which was in- 
serted into the container. In order to permit observation of phenomena 
occurring during discharge, a rotary Teflon valve was designed and installed 
in the test section base. When opened this valve allowed the liquid to be 
discharged through an unobstructed 1.6 mm wide slot. The heating section 
was also constructed primarily of Plexiglass. Aluminum plates were bonded 
to the inner vertical surfaces which had channels milled in them for water 
passage. These channels were connected to a water circulator and tempera- 
ture control bath. In operation, therefore, the exposed surfaces of the 
aluminum plates attained a uniform temperature very close to that of the 
circulating water, which was thermostatically controlled to + 0.1°C. 

Temperature measurements for reference purposes were made with four 
copper-aons tan tan thermocouples. Two of these were used to monitor the 
temperature in the bulk liquid near the bottom of the test section. The 
remaining thermocouples were placed in one of the aluminum plates at two 
different heights to permit confirmation of the temperature uniformity 
and measurement of the initial wall temperature transient. During a 
normal run, the aluminum side plates rapidly approached their steady state 
temperature after introduction of water frcm the circulator into the 



heating passages. The time constant for a 5°C wall temperature rise was 
approximately five seconds . The actual recorded variation with time was 
used as a boundary condition in the computer analysis. 


RESULTS 

Experimental data were obtained for values of T w (°°) - T q ranging up 
to 7°C. These data were extracted from infinite fringe interferograms 
such as the typical one shown in Fig. 2. A complete set of such photographs 
taken at intervals after the start of heating provides a detailed picture 
of the entire stratification process. One qu al itative observation of some 
Importance is the intense turbulent mixing action which occurs early in 
the process when the two streams of heated liquid moving along the free 
surface first merge at the center of the container. Laminar conditions 
were re-established fairly rapidly, but the temperature and flew fields 
were quite ooirplex during the next few seconds. 

A very interesting feature of the stratification process is the 
temperature reversal in the region near the wall. This is shewn clearly 
in Fig. 2. It can be seen that the temperature initially decreases with 
distance from the wall and then increases again. This is caused by two 
effects . 

(1) Since the Prandtl number is greater than one, the hydro- 
dynamic boundary layer is thicker than the thermal layer. 

Thus cool liquid from the bulk fluid moves upward along 
with the heated fluid within the thermal boundary layer. 

(2) The downward moving stratified liquid in the oentral region 
of the container is at a high temperature relative to that 
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of the upward moving liquid in the unheated portion of 
the hydrodynamic boundary layer. Thus the temperature 
reversal shewn results. 

Experimentally determined temperature profiles obtained from one set 
of interferograms at three different vertical positions in the container 
are shown in Fig. 3, together with calculated results obtained from the 
modified computer program. In order to illustrate the region of tempera- 
ture reversal clearly, the curves of Fig. 3 are drawn on an expanded 
scale. A test duration of two minutes was felt to be about the maximum 
over which accurate results could be expected due to heat losses from the 
container. At 123 seconds the largest error occurred in the lewer portion 
of the container. This error (based on the overall temperature difference 
between the wall and the bulk liquid) is only 23 percent. The agreement 
at other times and positions was much better. 

It will be noted that the experimental temperatures are lewer than 
the computed values in all cases. This is to be expected, since heat 
losses were present in the experimental system which were not accounted 
for in the analysis. In particular, energy was lost at the free surface 
of the liquid by conduction and evaporation, and at the glass end windows 
by conduction. The former effect was found to be negligible; an approxi- 
mate energy balance indicated that the energy loss to the glass windows 
accounted for most of the discrepancy between analysis and experiment . 

In addition to the quantitative measurements described, qualitative 
observations of the stratified fluid during discharge from the container 
were made. This was accomplished by taking movies of the interference 
pattern after suddenly opening the discharge valve. It was found that the 
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discharge process was basically laminar during the early stages with no 
significant mixing of the liquid. When approximately two thirds of the 
liquid had been discharged the process changed rapidly and an intense 
mixing occurred during the final phase. 


CONCLUSIONS 

The good agreement between analysis and experiment, particularly 
when the energy losses in the experimental system are taken into account, 
appears to verify the description of the stratification process provided 
by the work of Barakat and Clark. One effect present in an actual system 
which is not accounted for in their analysis is the turbulent mixing 
which occurs when the two streams of heated fluid first merge. The effect 
of this on the overall process at later times is obviously a function of 
the container size and the terrperature ratios involved, as well as the 
fluid properties. 
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Figure 1. Experimental Test Apparatus 
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Figure 2. Typical Interf erogram 
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USE OF MULTIPLE OUTLETS FOR REDUCING THE 
STRATIFICATION-INDUCED PROPELLANT SUPPLY TEMPERATURE RISE 
R. J. Schoenhals*, R. N. Nelson*, and E. R. F. Winter* 

Heat transfer to the fuel or oxidizer of 
a liquid fuel rocket produces thermal 
stratification, and this typically 
causes an increase in the propellant 
supply temperature during the latter 
portion of the engine burning period 
when the warmer fluid in the stratified 
layer is being discharged from the 
supply tank. This behavior is undesir- 
able because of the potential problem of 
cavitation in the supply pumps. In this 
paper a technique is described for 
reducing the stratification-induced 
supply temperature rise. Multiple out- 
lets are employed which feed into a 
mixing chamber below the supply tank. 

Experimental measurements are presented 
which illustrate the value of the 
proposed concept. 

INTRODUCTION 

During operation of liquid fuel rockets heat is unavoidably 
transferred to propellant supply tanks containing cryogenic 
fluids, both prior to and during flight. As the liquid in 
the boundary layer adjacent to the tank wall is heated, it 
rises due to buoyancy and is replaced by cooler liquid from 
the central core (Fig. 1)^. The upward -moving fluid moves 
radially inward as it approaches the free surface, and then 
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STRATIFIED LAYER 



LIQUID PROPELLANT 
SUPPLIED TO PUMP AT 
TEMPERATURE T p 

Fig. 1 Stratification in a Propellant Supply Tank 
Due to Heat Transfer 

decends slowly in the center. The temperature at the liquid 
surface increases as the process proceeds. Since warm fluid 
continues to rise along the wall, the heated fluid layer at 
the top becomes thicker. This layer is referred to as the 
stratified layer. The mechanism by which it is formed is 
called stratification. Detailed studies of this phenomenon 
are described in some of the references. 3 ^ 

A number of difficulties have been encountered in liquid 
propellant rocket technology which have been traced directly 
to the presence of stratification in supply tanks. Among 
these, one of the most serious involves the potential 
problem of cavitation in the supply pumps, 3,4,7 11 especial- 
ly near the end of the fuel-burning period when fluid in 
the stratified layer has descended to the bottom of the 
tank. This liquid is hotter than that discharged earlier in 
the burning period. Thus, its vapor pressure is considerably 
higher so that the onset of cavitation of this liquid, when 
it reaches the pumps, is far more probable. 3 ' 4,7 9 Once 
severe pump cavitation is established, the liquid flow may 
be choked off so that premature combustion burnout may occur 
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in the engine due to lack of propellant supply. 4,10,11 If 
this should happen when a substantial portion of the liquid 
still remains in the supply tanks, it is conceivable that 
the vehicle could fail to accomplish its mission. ' ' 

Or, on the other hand, even if mission accomplishment should 
result, it might be achieved at the expense of a considerable 
weight penalty if engine thrust is terminated due to 
cavitation with a large amount of unburned fuel and oxidizer 
still remaining in the supply tanks. A final possibility 
is that such cavitation can be prevented by increasing the 
supply pressure. This requires increased structural strength 
of the tank, however, so a weight penalty is incurred with 
this alternative also. 3,9,11,12 Nevertheless, this is the 
most commonly used procedure for avoiding cavitation. 

Several methods of eliminating the undesirable effects of 
stratification have been studied, and each of these possesses 
some disadvantages in terms of effectiveness, cost, weight, 
power requirements, complexity, and unreliability. The 
techniques suggested previously include the use of baffles 
and bubble lift pumps. 8,9 ' 12 

The present paper describes a different method for decreasing 
the stratification-induced propellant supply temperature 
rise during drainage from the supply tank. To the knowledge 
of the authors, this procedure has not been reported pre- 
viously. Experimental data are given to provide evidence of 
its future potential as a practical technique. A small 
scaled-down tank was used to obtain these measurements. 

Water was incorporated as the test liquid. The most favor- 
able measurements exhibit a rise in exit fluid temperature 
approximately 44% less than the corresponding rise which 
occurred during a conventional tank discharge process. This 
improvement indicates that a considerable decrease in 
maximum vapor pressure could be obtained for the liquid 
supplied to the pumps, thus enhancing the prevention of pump 
cavitation . 
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The proposed technique relies upon proper adjustment of the 
flow distribution for the liquid leaving the tank through 
multiple outlets. There is no requirement for additional 
pumps, power-driven mixers, actuation of valves during 
flight, or any other device containing moving mechanical parts 
and needing auxiliary power for its operation. Thus, the 
technique possesses certain inherent advantages of high 
reliability. It is clear, however, that careful engineering 
design, additional development, and increased complexity of 
construction would be necessary in order to produce favorable 
results with a full scale system. The following section 
contains a more detailed description of the method. 

INFLUENCE OF STRATIFICATION ON PROPELLANT SUPPLY 

Consider the stratification condition depicted in Fig. 1 as 
being representative of the situation just prior to engine 
ignition. The liquid can be considered to be initially at 
bulk temperature, T b , but heat transfer from the container 
wall has increased the temperature level of a portion of the 
liquid which has risen to the top. Here there is a stratified 
layer of liquid of thickness A which posseses the temperature 
distribution indicated in Fig. 1. Once ignition has occurred 
the liquid level decreases with time and, simultaneously, 
heat transfer to the liquid continues during the transient 
time period when the liquid is being discharged. 

If the total heat transfer occurring during the discharge 
period is small relative to that which has occurred prior 
to ignition, then the propellant supply temperature, T , 
varies with time in the same manner as T varies with height 
in Fig. 1. This result is illustrated by the solid line in 
Fig. 2. The propellant temperature first remains constant at 
the bulk temperature T b for a time t b before increasing. If 
all of the propellant is subsequently discharged from the 
tank, then the last portion of liquid leaving the container 
as it becomes empty (at time t g ) possesses temperature T g , 
the liquid surface temperature at the initial height H above 
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Fig. 2 Influence of Stratification on 
Propellant Supply Temperature 

the outlet, as shown in Fig. 1. If the heat transfer occur- 
ring during the burning period is appreciable, as it sometimes 
is when aerodynamic heating is involved, then the resulting 

T vs t curve lies somewhat above the one shown in Fig. 2 
P 

on the right hand side. However, this transient heating 
effect is omitted from the remaining discussion for purposes 
of clarity. The influence of the heat transfer during the 
discharge period can easily be incorporated as a later 
extension of the present study, since its occurrence does not 
disrupt nor in any way invalidate the general approach 
developed here. 

Now consider the situation in which severe pump cavitation 

occurs whenever the propellant supply temperature exceeds a 

cavitation value, T . If such a condition causes termination 

c 

of engine burning, then the total burning time is t instead 
of t as shown in Fig. 2. In this case a substantial portion 
of the initial propellant may be left in the tank without 
being consumed. This condition can be averted by increasing 
T to a higher value, T', by increasing the tank pressure. 

The higher pressure is sometimes achieved by employing a 
foreign gas (different from the vapor of the liquid 
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propellant) as a pressurizing agent. This usually requires 
a propellant tank with increased structural strength capable 
of sustaining the increased supply pressure, thus resulting 
in the weight penalty mentioned previously. 

A second technique involves the use of mixers or pumps which 

impart motion to the liquid in the supply tank to bring 

colder liquid from the lower portion of the tank into contact 

with the hotter liquid at the top to produce a mixture at an 

intermediate temperature level. If the liquid in the tank 

were completely mixed just prior to engine ignition, then the 

entire mass of liquid would be at a uniform temperature, 

T , the average temperature for the initial stratified 

temperature profile illustrated in Fig. 2. In this ideal 

case the propellant supply temperature Tp is constant 

throughout the burning period with T p = T & , and cavitation 

does not occur as long as the supply pressure is sufficient 

to make T > T . This desirable result is achieved only at 
c a 

the expense of the increased complexity associated with the 
incorporation of mixers, conventional pumps, bubble lift 
pumps, or other devices which can adequately prevent the 
hotter fluid from rising to the top and remaining there 
permanently. Baffles, which can be used to deflect the warm 
fluid rising near the tank wall toward the central core, 
offer a means of inducing some mixing without requiring the 
complexity of any mechanically moving parts nor relying upon 
the actuation of auxiliary valves. The use of baffles 
possesses several advantages over other techniques that have 
been considered, but this method suffers from the fact that 
the baffles are effective in preventing the stratification 
buildup only during the early portion of the transient 
stratification period. Further information on the techniques 
mentioned in this paragraph, including the advantages and 
disadvantages of each, can be obtained by referring to the 
appropriate references cited in the Introduction. The 
following section of the paper describes an alternate method. 
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Fig. 3 Propellant Tank with Two Outlet Tubes 

not previously reported, which forms the basis of the present 
study. 

USE OF MULTIPLE OUTLETS 

It is conventional practice to connect the propellant line 
or lines to the bottom of the propellant supply tank (Fig. 

1) . It is suggested here, however, that additional supply 
outlets can be used to remove the hotter fluid near the 
liquid surface continuously throughout the discharge process. 
Thus, the hottest fluid does not remain in the tank until 
the very end of the discharge process, nor is it necessary 
to employ a mixing device of some kind to reduce the undesir- 
able influence of the stratification on the propellant supply 
temperature . 

A schematic diagram of an initial design, incorporating only 
two outlets, is shown schematically in Fig. 3. With this 
arrangement. Valve A can be adjusted so that it constitutes 
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a sufficiently large flow restriction (relative to Valve B) 
so that an appropriate flow of hot fluid from the top is 
obtained. The hot and cold streams are brought together in 
the Mixer where an intermediate temperature is achieved for 
the single mixed fluid stream which flows to the pump. Two 
difficulties become immediately apparent at this point, and 
they complicate the detailed design of such a system: 

1. The flow of fluid from the upper outlet causes 
a variation in the temperature distribution 

in the stratified layer in the vicinity of the 
liquid surface. 

2. During the discharge process the upper outlet 
becomes completely inoperative as soon as the 
liquid level drops below the outlet. 

The first item appears to be an inescapable feature of the 

process which would need to be dealt with carefully in order 

to arrive at a final detailed design which derives the 

maximum possible benefit from the method. The second item 

suggests that it would be desirable if the mechanism could 

be made to operate throughout the discharge process. This 

can be accomplished by use of multiple outlets (more than 

two) located at various heights in the tank. Then, as the 

liquid level drops below the uppermost outlets, causing them 

to become inoperative, outlets at lower locations continue 

to operate and become effective as the hotter liquid is 

lowered into their vicinity. 

Although the use of several tubes having different lengths 
may give rise to desirable propellant supply temperature 
characteristics, the additional weight and the need to 
structurally support each of them within the supply tank 
appears to make this approach somewhat impractical. There- 
fore, in this study a compromise was sought which required 
only two outlet tubes, as shown in Fig. 3, but still pro- 
vided for multiple outlets at various heights inside the 
supply tank. This was accomplished by incorporating holes 
of different sizes along the length of the longer tube shown 
in Fig. 3 so that liquid was continuously removed at all 
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Fig. 4 Propellant Tank with Multiple Outlets and 

Two Outlet Tubes 

levels from the tank as illustrated in Fig. 4. The design 
of this kind of system involves the difficult question of 
the proper sizes of the holes in the tube, the number of 
holes to be used, and the distribution of the holes along 
the length of the tube. It was reasoned that the holes 
located along the upper portion of the tube should be the 
largest, with smaller holes further down. This arrangement 
should cause a sizeable flow of the hot fluid from the 
vicinity of the liquid surface very early in the discharge 
period, while producing a relatively larger restriction to 
the flow from the lower portion of the tank. Then, when 
the liquid surface descends and causes the uppermost holes 
to become inoperative, the holes further down become more 
dominant. Since there are holes allowing flow of the 
warmer liquid near the free surface at all levels along the 
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tube, this liquid is continuously removed and mixed with 
colder fluid from the lower portion of the tank. The fact 
that the flow restriction for the fluid near the surface 
increases when smaller holes are encountered as it decends 
is partially compensated for by the decreasing length of the 
flow path between the liquid surface and the Mixer. 

The ultimate goal would be the achievement of a propellant 
supply temperature which is essentially constant throughout 
the discharge period (as indicated by the dashed line in 
Fig. 2, for example). This may be possible with the proper 
size and distribution of the various holes along the 
discharge tube and, perhaps, even without the need for the 
second discharge line containing Valve A (see Fig. 4). This 
goal represents a problem of considerable magnitude and one 
which can probably be approached by means of a detailed 
analysis of the transient flow and thermal conditions in 
the tank during discharge of the liquid. Such an analytical 
study was not performed in the present research, however, 
sinae this investigation represented only a preliminary 
feasibility study of the technique in which complexity was 
avoided wherever possible. Therefore, a small experimental 
program was executed using the general arrangement shown in 
Fig. 4. Several different tubes with various hole configura 
tions were used, and Valves A and B were adjusted in an 
attempt to obtain a measured outlet temperature whose 
maximum value during the entire discharge process was as low 
as possible. These experimental results are described in 
the following section. 

DESCRIPTION OF EXPERIMENTS 

The experimental system contained a small tank 12 inches in 
diameter and 18 inches high fabricated from 1/8-inch thick 
brass. Water was used as the test fluid. The tank was 
filled with room temperature water to a depth of 16 inches 
prior to each test. Then, two 450 watt electrical strip 
heaters, firmly bolted to the outside of the tank, were 
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TEMPERATURE IN CENTRAL CORE , T (°F) 


Fig. 5 Typical Temperature Profile as Measured in the 
Central Core just Prior to Discharge 

turned on for a period of 15 minutes. The electrical power 
to these heaters was then terminated, and this was followed 
by a 3 minute waiting period to allow for the transient 
cooling of the heaters. Then a main control valve located 
below the Mixer was opened to initiate the discharge process. 
Thus, gravity driven flow was used with both the liquid 
surface and the outlet exposed to the atmosphere. This 
caused the flow velocity to decrease somewhat during the 
process. Any resulting deviation from the flow behavior 
associated with actual propulsion systems was deemed to be 
unimportant for purposes of this feasibility study. 

Some of the energy provided by the strip heaters was lost 
to the environment, some was utilized in heating up the 
tank itself, and the remainder was transferred to the water. 
However, as a result of the standardized 15 minute heating 
period and 3 minute waiting period, the temperature distri- 
bution in the liquid at the initiation of discharge was 
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essentially the same for all tests. This recurrent initial 
temperature distribution yielded a good basis for comparison 
of the various results obtained. Fig. 5 contains a typical 
plot of the temperature distribution in the central core, 
just prior to the initiation of discharge, as indicated by 
thermocouples located in the water at various depths beneath 
the surface. 

The lines leading to the Mixer were 5/16-inch diameter tub- 
ing, and the Mixer itself was merely a T-connection in this 
small apparatus. A thermocouple located downstream from the 
Mixer was used to measure the equivalent propellant supply 
temperature, T , of the mixed flow. The response of this 
thermocouple was displayed continuously by means of a strip- 
chart recorder. Several different supply tubes were mounted 
inside the tank and tested (Figs. 3 and 4). These tubes 
were fabricated from 1/2-inch diameter bakelite tubing. It 
was felt that this kind of tubing would prevent appreciable 
heat transfer from the hot fluid descending in the tube, and 
was therefore more appropriate than metal tubing for purposes 
of simplicity in this feasibility study. This would not 
necessarily be a requirement in a perfected design for an 
actual propulsion system, for indeed, it may be possible to 
achieve better performance by allowing the hot fluid to 
reject much of its excess energy to the central core before 
it reaches the Mixer. 

A sensitive pressure transducer was located at the bottom 
of the tank, and the signal from this transducer was also 
displayed on the strip-chart recorder. Since the water 
surface descended quite slowly during the discharge process, 
dynamic effects for the liquid flow in the tank were 
negligibly small. Thus, the variation of the pressure 
signal provided an instantaneous indication of the hydro- 
static pressure alone, which was proportional to the height 
of the water column in the tank. This system was calibrated 
so that a continuous plot of the liquid surface position as 
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Fig. 6 Typical Variation of Liquid Surface Position 

a function of time could be obtained during each test. 

Time plots of the liquid surface height had the 

general shape illustrated in Fig. 6 as a result of the 

decreasing flow rate associated with gravity driven flow. If 

the flow rate were constant, S would vary linearly with t 

and a plot of T p vs S would have exactly the same shape as 

a plot of T vs t. Since S did not vary linearly with t in 

these experiments (see Fig. 6), a plot of Tp vs S for a given 

test was slightly different from the corresponding plot of 

T p vs t, although the two plots were qualitatively very 

similar. It was found that plots of T p vs S were somewhat 

more useful than plots of T vs t. This was because signif- 

icant variations in T could be easily correlated with the 

P 

corresponding instantaneous relative position of the liquid 
surface with respect to the various openings in the bakelite 
discharge tube inside the tank. For this reason, the two 
strip chart recordings (T vs t and S vs t) were cross- 

r 

plotted to yield T vs S. 
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Fig. 7 Measured Outlet Temperature Response for a 
Conventional Discharge Process 

EXPERIMENTAL RESULTS 

Since the bulk water temperature, T fa , varied slightly among 
the tests, comparisons between tests are facilitated by the 
use of (Tp - T^) , rather than T p alone. Also, for the pur- 
pose of comparing the multiple outlet tests with a correspond- 
ing conventional test with drainage only from the bottom of 
the tank, this temperature difference is compared with a 
reference temperature difference, (T p Max - T^) , where T p Max 
is the maximum outlet temperature measured during the 
conventional discharge process. In all cases the minimum 
recorded value of T p was essentially equal to T^, the initial 
bulk temperature of the liquid in the tank. With the two 
temperature differences defined above, a convenient 
dimensionless basis for comparison is achieved with plots of 

(T - T. )/(T_. „ - T, ) as a function of S/H. Some typical 

pop Max o 

results are now presented in this manner. 

Fig. 7 shows the measured outlet temperature variation for a 
conventional discharge process, in the absence of an outlet 
tube inside the tank, with drainage occurring only from the 
bottom of the tank. As expected, the major features of the 
curve are qualitatively very similar to the plot of the 
initial temperature distribution, as indicated by the plot 
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Fig. 8 Comparison of Multiple Outlet Response with 
Conventional Response 


of T vs x in Fig. 5. There was a slight dropoff in T p near 
the end of the process. This dropoff was probably due to 
heat transfer from the hot surface to the air above it during 
the discharge period. The T p Max - value was found to be 
fairly close to the initial surface temperature, so the 
dimensionless temperature response contained in Fig. 7 may 
be loosely regarded as a plot of (T p - T b)/( T S “ T b^ * 

Figure 8 contains a comparison of the conventional response 
with a response obtained using a multiple outlet system. In 
this case there were only 3 openings in the upper discharge 
tube, one consisting of the open upper end of the tube 
located about 3 inches below the initial surface level, and 
the other two consisting of 3/16-inch diameter holes located 
about 6 1/2 inches above the bottom of the tank. The 
behavior indicated by the temperature response is easily 
interpreted in terms of the locations of these outlet 
openings. Early in the process there was considerable flow 
of hot fluid through the upper opening, and this caused the 
outlet temperature to rise. Then, the liquid surface 
dropped below this opening causing that flow to be terminated, 
and the outlet temperature rapidly dropped. When the hot 
fluid near the surface descended into the vicinity of the two 
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Fig. 9 Comparison of Multiple Outlet Response with 
Conventional Response 


lower openings the outlet temperature again rose temporarily, 
and then decreased again as the surface dropped below these 
openings. At the end of the process the outlet temperature 
still rose to the same maximum level as in the case of the 
conventional test, so evidently the flow of the hot fluid 
through the upper outlet tube was not sufficient to effect 
any improvement with regard to the maximum outlet temperature. 
Obviously the arrangement used did not produce any major 
improvement in this instance, but it will be shown that 
considerable improvement can be attained with a more favor- 
able distribution of the outlets. 

Figures 9 and 10 illustrate typical responses for upper out- 
let tubes containing many 1/8-inch diameter holes. Both 
of these cases exhibit substantial improvement over the 
conventional situation. Figure 9 shows a case in which the 
maximum outlet temperature rise was lowered by 41% relative to 
that of the conventional test. It appears from the response 
curve that the sizeable rise in T p near the end of the process 
probably could have been lowered if a substantial flow of 
hot fluid could have been sustained throughout the entire 
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Fig. 10 Comparison of Multiple Outlet Response with 
Conventional Response 

test. It is clear from the low values of T p recorded in the 
middle portion of the test that very little flow of the 
hotter liquid was obtained when the liquid surface was 
roughly halfway to the bottom. This can be seen by referring 
to the response curve of Fig. 9 in the range of S/H between 
0.45 and 0.65. 

Figure 10 shows the response curve for a similar situation, 
but with a somewhat more favorable arrangement of the 1/8- 
inch diameter holes in the upper discharge tube. Also, in 
this case the open end of the tube was beneath the liquid 
surface initially instead of being above it (see insets in 
Figs. 9 and 10) . A moderate flow of hot fluid was maintained 
throughout most of the process, but this did not eliminate 
the substantial temperature rise near the end that has 
characterized each of the responses previously discussed. 

The maximum rise indicated in Fig. 10 is slightly less than 
that exhibited in Fig. 9, however, which does represent a 
more desirable configuration for the outlet system. The 
maximum rise was 44% below that of the conventional test. 

This represents the most favorable performance achieved in 
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this study. It is evident from the response that the flow 
of the hotter fluid was quite small in the early portion of 
the process, but had a gradually increasing trend as the 
discharge progressed. 

It should be emphasized that the arrangements of the outlet 
holes used for obtaining the responses in Figs. 9 and 10 were 
established largely by trial and error, but with the 
application of some intuitive reasoning that developed in 
analyzing previous responses as indicated by the foregoing 
discussion. A total of about 25 tests were conducted, 
although only a few typical responses are presented here. 

It should also be mentioned that the total discharge times 
were not the same for all tests due to the different valve 
settings and different outlet hole configurations. Actually, 
only two different valve settings were used, and the drainage 
times for all tests were in the range from 6 to 12 minutes. 

The tests associated with Figs. 7 and 8 required approximately 
6 1/2 minutes, while those associated with Figs. 9 and 10 
required about 11 minutes. 

It appears that a desirable modification of the design 
associated with Fig. 10 should provide a greater flow of hot 
fluid in the early stages of the process and slightly greater 
flow in the intermediate stages. If such a design is achiev- 
ed, it may be possible to obtain a considerable attenuation 
of the sharp temperature rise which occurred near the end of 
each of the tests conducted in this study. Indeed, it may 
even be possible to achieve a nearly uniform value of T p 
throughout the entire process with a properly optimized 
design of the multiple outlet system. 

CONCLUSION 

It is concluded that the multiple outlet technique 
described in this paper is capable of effecting significant 
reductions in the stratification-induced propellant supply 
temperature rise. In this feasibility study the most favor- 
able experimental results exhibited a 44% decrease in this 
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temperature rise. These preliminary measurements indicate 
that improved outlet designs should make it possible to 
obtain further reductions in the maximum outlet temperature, 
but the refinement of such designs will require additional 
development effort. The proposed technique does not require 
auxiliary power, nor does it necessitate the use of any 
moving parts. Hence, the method is potentially capable of 
offering high reliability. If the technical advantages 
associated with this approach are sufficient to warrant 
further study on a practical basis, then analyses accounting 
for weight, cost and fabrication considerations would be 
desirable . 
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NOTATION 

H initial height of the liquid surface above the 

bottom of the supply tank 

S instantaneous position of the liquid surface 

relative to the initial liquid surface height 
t elapsed time since initiation of discharge 

t fa time at which propellant supply temperature begins 

to rise 

t time at which onset of cavitation occurs 

c 

t total time required to empty the supply tank 

T initial temperature in the central core 

T initial average temperature of the liquid in the 

supply tank 

T^ bulk temperature of the liquid in the supply tank 

T c propellant supply temperature at onset of 

cavitation 

T^ propellant supply temperature at onset of cavitation 
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under increased pressure 

Tp propellant supply temperature 

T d maximum value of the propellant supply temperature 

T g initial temperature at the liquid surface 

x distance from the bottom of the supply tank 

A initial thickness of the stratified layer 
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